Path dependent scaling of geometric phase near a quantum multi-critical point 
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We study the geometric phase of the ground state in a one-dimensional transverse XY spin chain 
in the vicinity of a quantum multi-critical point. We approach the multi-critical point along different 
paths and estimate the geometric phase by applying a rotation in all spins about z-axis by an angle t). 
Although the geometric phase itself vanishes at the multi-critical point, the derivative with respect 
to the anisotropy parameter of the model shows peaks at different points on the ferromagnetic side 
close to it where the energy gap is a local minimum; we call these points 'quasi-critical'. The value 
of the derivative at any quasi-critical point scales with the system size in a power-law fashion with 
the exponent varying continuously with the parameter a that defines a path, upto a critical value 
a = Qc = 2. For a > Qc, or on the paramagnetic side no such peak is observed. Numerically 
obtained results are in perfect agreement with analytical predictions. 
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I. INTRODUCTION 

Quantum Phase Transitions (QPT) in quantum many 
body systems are being studied extensively in recent 
yearsir— . The possibility of experimental studies on ul- 
tracold atoms trapped in optical lattices, which can for 
e.g., undergo a mott insulator to a superfluid transition, 
have opened new avenues to investigate quantum phase 
transitions^!^. QPTs occur at absolute zero tempera- 
ture and are associated with a fundamental change in the 
ground state of the system^. A quantum critical point 
is characterized by a diverging length scale as well as a 
diverging time scale, namely the relaxation time of the 
quantum system. This characteristic time scale is the in- 
verse of the minimum energy gap of the underlying quan- 
tum Hamiltonian which vanishes at the quantum criti- 
cal point. Recently, QPTs have also been studied from 
the viewpoint of quantum information theory. Quanti- 
ties like concurrence^, entanglement entropy^i^, fidelity 
susceptibilitj«2i"— and geometric phasesi^r^ have been 
studied in the vicinity of quantum critical points, and arc 
found to capture the singularities associated with them. 
A close relation between geometric phase(GP) and mag- 
netic susceptibility has also been established for trans- 
verse XY chain for QPTs driven by transverse magnetic 
field^^. It is worth mentioning that the ground state GP 
in a Heisenberg XY Model has been experimentally ob- 
served very recently using NMR interferometry^l. 

In the present work, we investigate the scaling behav- 
ior of GP— near a quantum multi-critical point (MCP) 
generated by the application of rotation on each spin of a 
quantum spin chain around z-axis. Berry showed in that 
in addition to the usual dynamic phase, a phase purely 
of geometric nature is accumulated on the wave func- 
tion of a quantum system under an adiabatic and cyclic 
change of the Hamiltoniar^. Quantum critical points 
are regions of vanishing energy gaps, and consequently 
are accompanied by non-analyticities in various observ- 



ables of the system. In a quantum critical system, the GP 
of the ground state of the system, which depends on the 
energy gap, can capture the associated singularitiesi^"— . 
In fact, GP can also be related to the imaginary part of 
a general geometric tensor, whose real part on the other 
hand gives the fidelity susceptibility which is the rate of 
change of the ground state of the Hamiltonian following 
an infinitesimal change in its parametersi2iiii21i^. 

We exploit the integrability of spin-1/2 XY chain^^Ti^^ 
to investigate the nature of the geometric phase close to 
the MCP. Recently, it has been shown that a GP dif- 
ference between the ground and the first excited state 
exists in an isotropic XY chain (i.e., the XX chain) if 
and only if the closed evolution path circulates a region 
of criticalityii^. In this paper, we use a dynamical scheme 
so that the MCP is approached along various paths char- 
acterized by a parameter a. The non-contractible GP of 
the ground state is studied following the behavior of GP 
and its derivative close to the Ising transition point of an 
anisotropic XY chair^i^. In the case of a MCP, we find 
that the derivative of the GP with the anisotropy param- 
eter scales with the chain length in a non-trivial fashion 
with an exponent depending on the parameter a. More- 
over, it does not peak right at the MCP, rather it peaks at 
points close to the MCP on the ferromagnetic side (which 
are so-called 'quasi-critical points— i^), where the energy 
gap is a local minima. However, beyond a limiting value 
oi a ~ ac no such peak is observed and GP as well as its 
derivative at the MCP is zero for all values of a. 



II. THE MODEL AND THE GEOMETRIC 
PHASE. 

The model we consider is a one-dimensional spin- 
1/2 XY model in a transverse field with near- 
est neighbor ferromagnetic interactions given by the 
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Hamiltonian29|30^^ 

^ = -^ E [(i+7KX+i + (i-7KX+i + /^<],(i) 

n 

where cr's are Pauh spin matrices satisfying the usual 
commutation relations. The parameter h is the mag- 
netic field applied in the z— direction and 7 measures the 
anisotropy in the in-plane interactions. The phase dia- 
gram of the system, plotted in /i — 7 plane is shown in 
Fig. ^ where the vertical bold lines at /i = ±1 repre- 
sent quantum transitions from the ferromagnetic to the 
paramagnetic phases that belong to the transverse Ising 
universality class and hence called the 'Ising transitions'. 
The horizontal bold line at 7 = with the transverse field 
lying between 1 and —1 represents transitions between 
ferromagnetic ordered phases with ordering in x and y di- 
rections, respectively^!^. The points A and B (/i = ±1 
at 7 = 0, respectively) where Ising and anisotropic lines 
meet are the multi-critical points^^ which are our sub- 
ject of interest. Analyzing the energy spectrum of the 
Hamiltonian ([Tp^, it can be easily shown that the en- 
ergy gap scales with the momentum k as (= k) at any 
Ising critical point {h = ±1,7 7^ 0) so that the dynami- 
cal exponent z = I. On the other hand, for the critical 
mode the gap scales as \h ~ ll^^i— \h ~ 1|) yielding the 
correlation length exponent = 1. Similarly, one finds 
for MCP that z = 2 and v = 1/3^; i/z = 1 in either 
cases. 
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FIG. 1: The phase diagram of a one-dimensional XY model 
in a transverse field. The vertical bold lines {h = ±1) denote 
Ising transitions from the ferromagnetic to the paramagnetic 
phase. The horizontal bold line stands for the anisotropic 
phase transition from a ferromagnetic phase with magnetic 
ordering in x-direction to a ferromagnetic phase with order- 
ing in y-direction.The multi-critical points are at A (ft = 1, 
7 = 0) and B (/i = —1, 7 = 0). We show different paths of 
approaching the MCP A corresponding to different values of 
a; path I (path II) is for a = 1 (a = 2). We show that for 
a < 2, we observe quasi-critical peaks and there is a continu- 
ously varying effective scaling exponent for dfigldrf with the 
chain length, while for a > 2, no peak is observed. 



We propose a dynamical scheme that enables us to 
approach the MCP along the path defined by^ 

Kl) = 1 - l7rsff"(7), (2) 

so that the system hits the MCP A at ft- = 1 , 7 = 
(see Fig. ([1])). Close to a MCP, the path of approach 
plays an important role as also to be shown below in the 
scaling of any dynamical response; the dynamical scheme 
(PI in which the path can be changed through tuning the 
parameter a therefore facilitates the study on the path- 
dependent scaling of CP close to the MCP. 

In order to investigate the GP in this system, we in- 
troduce a new family of Hamiltonians that can be de- 
scribed by applying a rotation of 77 around the z di- 
rection to each spini^iiS, i.e., if,, = grjHgjj with g,^ = 
n„ exp(i77cr^/2). For the family of Hamiltonians thus 
generated, the energy spectrum remains the same leav- 
ing the critical behavior unaltered. The Hamiltonian 
Hiji) can be diagonalizcd by using the Jordan- Wigner 
transformation which transforms the spin operators into 
fermionic operators a„ and at via the relations^Si^ 
= iY\i<n'^i)'^n- Employing a Fourier transformation 
dk — [aj cxp{—i2Trjk/N)] /y/N followed by Bogoli- 
ubov transformation, one can recast Hamiltonian ([T]) to 
the diagonal form H = J2k^k{cl.Ck — 1), where A^ = 
\/{h + cos(27rfc/A^))2 + 72 sm^(2TTk/N), Ck = dk cos ^ - 
d^j.e^'"' sin and the angle 6k is given by cos 9k = 

{cos^ + h)/Ak. 

The ground state is a tensor product of states, each 
lying in the two-dimensional Hilbcrt space spanned by 
|0)fe|0)_fe and |l)fc|l)_fc, where |0)a; and |l)fe arc the vac- 
uum and the excited state for k-th mode |l)A;|l)_fc = 
dkd''_f.\0)k\0)-k- Under the unitary transformation 
(i.e. rotation about z— axis), the ground state of the 
transformed Hamiltonian picks up an additional phase 
factor and is given by \g) = Hfclcos ^|0)fe|0)_fe - 
le^") sin ^|l)fc|l)_fc). The GP of the ground state, accu- 
mulated by varying the angle r/ from to tt, is described 
by2^ /3g = -i/N {g\drj\g)d-q , and is given by 

^ = ^T.^^-^^^^k). (3) 

k 

This equation in the thermodynamic limit {N — >■ 00), 
where the critical properties are studied , is given by 

Pg^ j {1- COS 0^)d<l>, (4) 
Jo 

where the summation 1/Nj2k replaced by the integral 
I/tt with (j) = 2TTk/N, so that 

cos 6*0 ~ (cos0 + /i)/A0 and 

A^ = J{h + cos 0)2 + 72 sin^ (j). (5) 
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Carollo and Pachos^ studied the behavior of the GP 
close to the anisotropic critical point and showed that 
a non-contractible geometric phase difference between 
the ground state and the first excited state exists when 
the Haniiltonian encounters a critical point while pass- 
ing through an adiabatic cycle. In a subsequent work, 
Zhiii^ showed that the derivative djSg/dh as h is varied 
shows a peak right at the Ising critical point and di- 
verges logarithmically with the chain length. From the 
scaling relations d/Sg/dh ~ KilniV -I- Ci at the critical 
point h ^ he = I, and d/Sg/dh ~ K2lD.\h — hd + C2 for 
an infinite system (where Ci and C2 are non-universal 
constants) , the critical exponent v can be obtained from 
the relation v ~ H2/ ni = 1- 



III. RESULTS 

We study the behavior of GP in the vicinity of the 
MCP 'A' (see Fig. ([T])) using the dynamical path given 
in Eq. Let us first consider the case of a = 1, i.e., 
a linear path approaching the MCP and estimate the 
geometric phase and its derivative using numerical tech- 
niques. We observe a series of peaks in the derivative 
dPg/d'y when plotted against 7 close to the MCP on the 
ferromagnetic side (7 > 0, \h\ < 1 for the path ^) of it 
as shown in Fig. ([2]). A similar behavior is observed for 
all values of a < 2 (see Fig. ^). This is in contrast to 
the behavior near an Ising critical point where there is 
only one peak right at the quantum critical point^^ (see 
Fig. dU). No such peak is observed on the paramagnetic 
side (7 < 0, \h\ > 1) and the GP itself becomes trivial 
at the MCP as seen from Eq. ([3]). On the other hand, 
for a > 2 no peaks are observed close to the MCP (see 
Fig. Q); and the case a = 2 shows a limiting behavior to 
be discussed later (see Fig. ([5])). The appearance of the 
series of peaks in d^g/dj can be attributed to the exis- 
tence of quasi-critical points where the energy spectrum 
attains a local minima, to be explained below. 

The scaling behavior of the derivative d/Sg/d-f with the 
system size for a generic a can be derived in the following 
way. From Eq. ^ we get the spectrum as 



A, 



\J [cos + (1 - |7|°s5n(7))]2 + 72 sin^ cj), 



so that using Eq. (|4]) for positive 7, dfUg/d-f can be ex- 
pressed as 



d^ 
c?7 




(7) 



where 

A = {cos(/)+(l-7")}{[cos(/)+(l-7")]a7""i-7sin2 cj,}. 
The first term of the integrand dominates over the sec- 
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FIG. 2; (Color Online) The numerically observed variation of 
d/Sg/d-y with respect to 7 for a = 1 for different chain lengths, 
iV = 2000 (blue solid line), = 1800 (red dotted line) and 
N = 1500 (green dashed line). In all cases, a series of peaks 
is observed on the ferromagnetic side (7 > 0, \h\ < 1) along 
the path of approach. With increasing A'' the peak height 
increases, the peaks shift closer to MCP, and the distance 
between subsequent peaks decreases. 
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FIG. 3: A similar behavior as in Fig. ([2]) for d/Sg/d'y as a 
function of 7 is observed for a = 1.5. 



ond term which can be dropped for further calculations. 
For ~ TT, the minima in A^ occurs at points where 
(/)^ (shown below) so that 
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(6) dl3^^ r 
d'-f 



;^)(i+2) 



with (p ^ l/N. Hence we get the scaling relation 



d^ 
c?7 



oc A^- 



(8) 



(9) 



The scaling relation ([9]) shows that dfig/dj diverges 
with the system size for a < 2 while for a = 2 a satu- 
ration is expected. The scaling exponent varies contin- 
uously with the path, i.e., with the parameter a upto 
the limiting value a = 2. Fig. ([6]) shows the scaling of 
d/3g/d'y with A'^ for a = 1 and 1.5 confirming the ana- 
lytical scaling ^ while no peak is expected for a > 2 
as numerically observed in Fig. With increasing a, 
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FIG. 4: Numerically no peak is observed in dj3g/d'y for a = 2.5 
which is due to the absence of quasi-critical points along the 
path. However, a peak is observed corresponding to the Ising 
critical point at 7 = 1.319, h = —1. 



FIG. 6: Scaling behavior of d/Bg/d"/ with A'^ is shown for a = 1 
and a — 1.5 (inset). The numerical results perfectly match 
with the scaling relation ([9|. 
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FIG. 5: The numerical plot of dj3g/d'y with 7 for the marginal 
case a = 2 indicating a peak close to the MCP (left inset) 
with the peak height saturating to a constant value with A'' 
(right inset) as expected from the analytical scaling relation 
The Ising peak is observed at 7 = 1.414 where h — —1. 



the divergences seen in Eq. (jH]) is expected to be slower 
explaining observations in Figs. ^ and 

Let us now focus on the energy spectrum ([HI); the 
energy gap is not minimum at the MCP for a < 2, 
rather it is minimum at the 'quasi-critical' points de- 
termined by the condition cost/) + (1 — \^\°'sgn{'y)) = 
or |7|"sgn(7) = (jp if we expand </> ~ tt and rescale 
(/) — > TT — 0. The derivative of GP shows a peak when- 
ever the energy gap is minimum, i.e., the spectrum hits a 
'quasi-critical' point resulting in a series of peaks on the 
ferromagnetic side of MCP as numerically observed in 
Fig. (|2|) and Fig. The condition |7|"s.gn(7) = 0^ docs 
also imply that for increasing chain length, the quasi- 
critical points (and hence the peaks in d^g/d'y) become 
closely spaced. Note that strictly for N ^ 00, the peaks 
collapse to a single peak, however for any finite chain, 
there will be multiple peaks. In fact even in the iV 00 
hmit, the quasicritical exponents (not the exponents as- 
sociated with the MCP) will appear in the scaling rela- 



tion ((9)) as also observed in the scaling behavior of defect 
density following a slow quench^ and that of the fidelity 
susceptibilitjs^. The analysis of Eq. (jH]) also shows that 
no such quasi-critical point exists on the paramagnetic 
side and hence no peak in djig/d'y is expected. Similarly, 
for a > 2 the minimum in energy gap occurs right at 
the MCP where the GP and its derivative vanishes and 
hence no peak appears as shown in Fig. For a = 2 
only one peak of d/3g/dj is observed close to the MCP 
though the height of the peak saturates for large N to 
a constant value as the exponent {4/ a — 2) in Eq. ^ 
vanishes in the limit at a — > 2 (Fig. ([5])). Investigating 
figures ([4]) and ([5]) closely, we observed peaks at the Ising 
critical points h = —1; this peak appears for all values 
of a and scales logarithmically with N yielding = 1 as 
reported earlier—. 

From Eqs. one finds that for a finite chain 

with isotropic interaction (7 = 0), 6^ = or tt so 
that /3g = or 2tt for h > 1 whereas for h < 1, 
fig = 2tt — 2 arccos(/i)^'2i. Additionally, in the thermo- 
dynamic limit, for infinitesimally small 7, we can always 
find a solution (po such that cos 00 4- (1 — \'y\"sgn{'y)) = 

but A0„ = 7v/l - (1 - |7|"S5'^(7))^ ^ (see Eq. jH)). 
This leads to 6^g = tt/2 and and a non-trivial geometric 
phase /3g = TT for 7 — > 0. 

The relation between the geometric phase and the 
transverse magnetization (AI^) at zero temperature given 
by the relation f3g = n + ttMz, has been established 
close to the Ising transitions^. Close to the MCP 'A', if 
one defines an a-dependent magnetization M{a) through 
the derivative dh-^/dh for a given path (see Eq. 
it can be shown that M{a) = + dA^/dj{dj/dh){= 
1/A^[{h + cosfc) — (2/a)7^~" sin^ fc]), where the trans- 
verse magnetization is related to the geometric phase 
as discussed in the reference^. One can therefore con- 
clude that close to a MCP the derivative of the geometric 
phase along a path can not be related to the transverse 
magnetic susceptibility directly. 
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IV. CONCLUSION 

In this work, the scaUng properties of GP have been 
studied close to a quantum MCP of a spin-1/2 transverse 
XY spin chain. We show that the scahng of the deriva- 
tive of the GP depends on the path of approaching the 
MCP. It shows a peak whenever the system hits a 'quasi- 
critical point' leading to a series of peaks on the ferromag- 
netic side of the MCP. The peak diverges with the system 
size with an exponent that varies continuously with path 
up to a limiting path given by a = ac = 2 where we ob- 
serve a saturation with A^. On the paramagnetic side and 
also for paths beyond the limiting path (i.e., for a > 2), 
the system does not encounter any quasi-critical point 
causing the peaks in the derivative to disappear. 

We note that in a recent study, a similar behavior 
of the fidelity susceptibility, namely, occurrence of se- 
ries of peaks on the ferromagnetic side of the MCP has 
been reported^. The fidelity susceptibility has also been 
found to diverge with the system size with an exponent 
that continuously varies with the path up to a = ac = 2. 
Our study therefore points to a deep connection between 
the scaling of the fidelity susceptibility and GP close to 
quantum critical as well as multi-critical points, as was 
predicted in the study of quantum geometric tensor^iSiii. 
Questions may remain on the nature of the quasi-critical 



points, especially whether they are generic or specific to 
this particular spin chain. These points may appear in 
the present model due to the existence of the line of 
anisotropic critical points with continuously varying or- 
dering wave vectors. However, recent studies including 
the present work reveal that in all dynamical responses 
near a quantum MCP, e.g., defect generation following 
slow and rapid qucnches^i^, fidelity susceptibilit}*^ as 
well as in the scaling of the GP, these quasicritical points 
play a dominant role depending on the path of approach. 

Clearly, the scaling presented in Eq. © and the limit- 
ing value a ^ ac{= 2) are established close to a MCP in 
a transverse XY chain. For a generic MCP, the scaling 
law and the value of ac or the position of the peaks are 
expected to be different as they depend on the spectrum; 
however near a generic MCP if quasicritical points ex- 
ist, the associated exponents shall dictate the equivalent 
scaling relations and the limiting path, and the derivative 
of the GP will show peaks at the quasicritical points. 
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